Atom-wall dispersive forces: a microscopic 

approach 

F. Cornu 

Laboratoire de Physique Theorique, UMR 8627 du CNRS 
Universite Paris-Sud, Bat. 210 
F-91405 Orsay, France 



Ph. A. Martin 
Institute of Theoretical Physics 
Swiss Federal Institute for Technology Lausanne 
CH-1015, Lausanne EPFL, Switzerland 

November 5, 2009 



Abstract 

We present a study of atom-wall interactions in non-relativistic quan- 
tum electrodynamics by functional integral methods. The Feynman-Kac 
path integral representation is generalized to the case when the particle 
interacts with a radiation field, providing an additional effective potential 
that contains all the interactions induced by the field. We show how one 
can retrieve the standard van der Waals, Casimir-Polder and classical Lif- 
shiftz forces in this formalism for an atom in its ground state. Moreover, 
when electrostatic interactions are screened in the medium, we find low 
temperature corrections that are not included in the Lifshitz theory of 
fluctuating forces and are opposite to them. 
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1 Introduction 



1.1 Issue a stake 

Atom-wall interaction plays an important role in several physical, chemical and 
biological situations [T| [2] and has been the subject of active measurement in- 
vestigations during the last years (see [3] for a review and |4J for references to 
more recent experimental results). The basic model consists of a quantum atom 
in its ground state at a distance X from a perfect metallic wall : the atom inter- 
acts both with its wall mirror image and a source-free radiation field. In 1932 
Lennard- Jones considered the case where interactions are purely electrostatic 
and predicted a l/X 3 atom- wall interaction [5]; the latter is similar to the van 
der Waals interaction between two atoms in their groundstates, first estimated 
by London in the framework of quantum mechanics and classical electrodynam- 
ics [6]. 

In 1948 Casimir and Polder adressed the case where the electromagnetic 
radiation field is quantized and they calculated the atom-wall interaction at 
the second order of the ground-state perturbation theory [7]. For X small (a 
few nanometers) they showed that the atom-wall potential is the electrostatic 
dipolar l/X 3 van der Waals-like interaction, with an amplitude given by the 
ground-state position fluctuations inside the isolated atom. At larger distance 
(a few micrometers) they found a cross-over to a longer-ranged l/X 4 potential, 
with an amplitude determined by the vacuum fluctuations of the source-free 
radiation field confined by the metallic wall. When calculations are performed 
in the Coulomb gauge of electrodymanics, the interpretation in terms of the 
typical atomic-line wavelength A a t is as follows. At distances X <C A a t the 
Coulomb image interaction dominates the atom-radiation interaction (of relative 
order A/A a t). However at distances X ^> A a t the atom-radiation interaction 
becomes of the same order as the Coulomb interaction. In fact its leading 
term becomes opposite to the Coulomb image interaction, while its subleading 
term (of relative order A a t/A) does not involve the instantaneous interactions 
artificially introduced by the Coulomb gauge : the 1/A 4 Casimir-Polder tail 
proves to be a retardation effect. 

For an atom prepared in an excited state and which interacts with the pho- 
ton field in its vacuum state, the atom-wall interaction has been investigated 
at the second order in the atom-field coupling in Ref.[8j. The authors used 
the microscopic approach of Ref.|9l ITU] (see also [11]) originally introduced to 
disentangle the two phenomena that come out in the source-field picture: on 
one hand the "self-reaction" effect, due to the interaction between the electron 
and its own field, on the other hand, the "vacuum fluctuation" effect, due to 
the interaction between the electron and the source-free part of the quantized 
electromagnetic field. The short-distance interaction is again dominated by the 
van der Waals-like interaction, which emerges as a purely self-reaction contri- 
bution. At large distances the leading net atom-wall interaction becomes an 
oscillating l/X decay : this tail is similar to the interaction between a classical 
oscillating electric dipole and its own reflected far field. However it arises only 
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from lower-lying states, as in the spontaneous-emission rate, because it results 
from the combination of self-reaction and vaccum-fluctuation contributions. For 
an atom in its ground-state, which cannot spontaneously radiate, the oscillating 
X/X tail proves to vanish and the atom- wall interaction is reduced to the sub- 
leading contribution, namely the Casimir-Polder potential which results from 
the change in the spectrum of vacuum fluctuations due to the reflection of the 
source-free radiation field on the metallic boundary. 

In the case of an atom in equilibrium with a thermalized photon field, the 
atom-wall interaction has been studied for a simplified model of the atom in 
Ref.|12] by using the microscopic approach of Ref.|9j [10] again. The atom in 
the model contains only two atomic levels, and the limitations of this simplica- 
tion have been pointed out in [13]. The energy shift of the excited state has an 
oscillating 1/X tail, but, after thermal equilibrium average, this contribution is 
exponentially damped at low temperature. Moreover, at distances larger than 
the length scale A p h beyond which thermal fluctuation contributions dominate 
over vacuum ones, the energy-shift X-dependence is changed into the 1/X 3 
classical Lifshitz potential [13] which arises only at non-zero temperature. In- 
deed, nowadays the semi-macroscopic Lifshitz theory |14l Il5] (see also [2]) is 
the commonly adopted theoretical framework for calculating forces induced by 
quantum and thermal fluctuations. A detailed application of Lifshitz theory to 
the atom-wall interaction can be found in the recent paper [3] (see references 
therein). 

In the present work, we investigate the low-temperature atom-wall interac- 
tion for an atom in a dilute gas in equilibrium with the photon field when the 
internal structure of the atom is described from the full Hamiltonian of quantum 
electrodynamics with nonrelativistic matter. 

We revisit the original Casimir-Polder model with a method based on func- 
tional integration which is in principle not perturbative and not limited to zero 
temperature. The model consists of an Hydrogen atom with spinless electron 
of charge e and mass m, and an infinitly massive proton. For a single atom 
in thermal equilibrium with the electromagnetic field at inverse temperature 
P = (fcs? 1 ) -1 , the atom-wall potential Q(X, (3) is defined as the excess free en- 
ergy of the atom immersed in the photon field at a distance X from the wall. The 
calculation of Q(X, j3) involves the trace over the photonic degrees of freedom 
of the quantum Gibbs factor. The main tool is the functional representation of 
the thermal equilibrium weight of a quantum particle submitted to an external 
static potential by means of the Feynman-Kac-Ito formula |16l W7\ . When the 
external field fluctuates, the integration over the fluctuations leads to a gener- 
alized Feynman-Kac-Ito formula with an additional effective potential. Using 
a bosonic functional integral representation of the thermalized photon field, we 
establish such a generalized Feynman-Kac-Ito formula for the electron in in- 
teraction with the quantum electromagnetic field in some region with general 
boundary conditions (see (26 1 and (28 1). In an homogeneous system where the 
radiation field extends over the whole space, we recover the results of [18] for a 
classical field and [19 for a quantum field. 
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1.2 Finite-temperature effects 



As soon as the temperature is different fron zero, numerous effects come into 
play and we have to investigate them separately. Among them we distinguish 

1. The spectral broadening and level shifts of the atom due to radiative 
interactions with the photons. 

2. The thermal excitation of the atom and its possible ionization. 

All these effects are embedded in the effective potential Q(X, (3). In this paper 
we determine the dominant term of Q(X, (3) for X large compared to various 
microscopic lengths in the system, while keeping only lowest order effects in the 
fine structure constant (%. This will amount to neglect spectral broadening and 
level displacements in final calculations and deal with the usual energy levels 
Ei, i = 0, 1, ... of the bare Hydrogen atom. Moreover it also corresponds to 
work in the so-called dipolar approximation (as in the original Casimir-Polder 
paper), where the radiation-field spatial fluctuations inside the atome are ne- 
glected. The study of relativistic diamagnetic terms, that arise beyond the 
dipolar approximation, is postponed to a forthcoming paper |20] . 

Concerning the second point, one must observe that the notion of atomic 
bound state makes sense only if the thermal energy is much less than the ion- 
ization threshold |i?o|, namely fc^T <C l-Ebl- We shall therefore consider the 
low-temperature regime characterized by 

k B T « Ei - £ (1) 

when the thermal energy is insufficient on the average to excite the atom from 
its ground state. More precisely, we shall neglect exponentially vanishing tem- 
perature corrections 0{e~^ El ~ E ° > ) and retain only those that are algebraically 
growing with T. 

At the microscopic level, even if the low-temperature condition is sat- 
isfied, there is always a non-vanishing probability to populate Rydberg states 
and to ionize the atom. As a consequence a summation on the atomic thermal 
weights is divergent. 

In reality envisaging a single thermalized atom does not make physical sense 
: one has to consider a dilute gas of atoms at positive density. In addition to 
thermal photons, at any non-zero temperature the atoms in the gas are also in 
equilibrium with an ionized fraction of electrons and nuclei that provide collec- 
tive screening mechanisms on the length scale A scree n- For a purely Coulombic 
quantum gaz at low density and low temperature, it has been shown that col- 
lective screening has two main consequences : 

(i) a natural regularization of atomic traces that makes them finite, 

(ii) the introduction of an appropriate screened Coulomb potential. 

These results come from the elaborate discussion of the subtle interplay be- 
tween quantum mechanical binding, ionization and screening presented in [5T] 
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as well as from the analytical investigation of the effective interactions in a 
partially-ionized dilute gas at low temperature in |22j or the effective charge- 
wall interaction in an ionized dilute gas [23]. 

Although we presume that the same studies extend to the present electro- 
magnetic system with the same conclusions, we shall not enter here in the full 
many-body problem in order to keep the presentation as simple as possible. For 
the finiteness of traces (point (i)), we rather introduce a formal spatial cut-off 
when calculating the traces of atomic observables by limiting the spatial integra- 
tion over the electron position r to a ball of radius i?o centered around the proton 
position R, f dr — > J rfr| r -R|<fl - This cut-off, of the order of R = p -1 / 3 (p 
the density of the screening gas), can be interpreted as delimiting the effective 
available configurational space for an atom when the density is different from 
zero. Equivalently, when working in the energy representation, we shall cut the 
level sums at some maximal energy J? max . AVhen X Ascreen i we calculate 
3>(X, 0) in the low-temperature regime (JlJ using the bare Coulomb potential. 
It turns out that in this regime, results are independent of the cut-off Ro which 
can eventually be removed. When X 3> A scree n we procceed in the same way, 
but take collective screening into account by replacing the bare Coulomb po- 
tential between the atom and its image by the screened potential mentioned in 
point (ii) above. This screened potential, extensively studied in |2l], does not 
follow the classical exponential Debye law, but it has an algebraically decaying 
tail. This tail, generated by the intrinsic quantum fluctuations of the particles 
|25] . is dominating when X ^> A scree n- It is of dipolar type and therefore should 
not be omitted in the calculation of the atom- wall interaction [23 . 

The X-behaviour of the atom-wall interaction is also determined by the 
hierarchy of the length scales in the microscopic model. First, since the electron 
is non relativistic, one has to disregard high-energy photons that could generate 
high-energy processes such as pair creation, namely we consider only photons 
with hck <C me 2 (k is the photon wave number and c the speed of light). This 
introduces an ultraviolet wavenumber cut-off k cut = (A cu t) 1 = mc/h on the 
field modes. Next, the length scales associated to the atomic properties are 
the Bohr radius a B = ti 2 /me 2 and the atomic wavelength A a t = hc/{E\ — Eq) 
corresponding to the photonic transition between the ground state and the first 
excited atomic level. Finally A p h = (3hc defines the thermalization length of the 
photon. Since E\ — Eq = 3|J3o|/4 and \E \ — e 2 /2a B , one has A cut /a B = ofo 
and a B /A at = (3/8)af 8 where otf a = e 2 /he is the fine structure constant, so that 

A cut < a B < A at (2) 

Since A a t/A pn = ksT /{E\ — Eq), the condition ([lj implies also that A at -C A pn 
and with Q the full hierarchy of length scales reads 

A C ut < a B < A a t < Aph (3) 

However it will be always assumed that the wall separation is larger than the 
atomic size, X ^> a B (the model does not make sense for X ~ a B since the wall 
is described at the macroscopic level). 
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1.3 Scheme of the paper and statment of results 

The paper is organized as follows. First, in section 2 the generalized Feynman- 
Kac-Ito formula with the additional effective potential that embodies the inter- 
actions induced by the quantized field is derived for general boundary conditions. 
The quantum photonic oscillators are represented by the Gaussian oscillator 
stochastic processes. In this representation, the electron appears as submited to 
an external classical random field and the standard Feynman-Kac-Ito formula 
applies. Then the field degrees of freedom can be integrated out by means of 
a simple Gaussian integral. The procedure is well known and can be found in 
[17] and [26] page 187. 

In section 3, we specify the above formula to the atom-wall model by in- 
troducing the mirror charges and fixing the field boundary conditions at the 
metallic wall. The expression of the atom-wall potential 3>(X, 0) in terms of 



path integral is given in (58 1. At this point, the large-distance analysis is not 



perturbative with respect to the dimensionless coupling constant af s . For in- 



stance, formulae (60) and (61 1 contain several higher-order effects in at a such 
as diamagnetic polarization terms and thermal displacement and broadening of 
spectral lines of the atom in the photon field. 

In section 4 we consider a density regime where A screen is larger than all 
length scales in the hierarchy ^ and fix the temperature in the low-temperature 
regime (JTJ where exponentially-small temperature corrections are negligible 
compared to the algebraic ones. We show how our formalism allows one to 



recover the van der Waals (95 I, Casimir-Polder (99 1 and classical Lifshitz (1031 
potentials in the respective ranges X <C A at , A at <C X <C A p h and A p h <§C X <C 
Ascreen- For this aim we neglect the above mentioned finer relativistic effects by 
making the dipolar approximation and by switching off the atom-field coupling 
in the bulk. Then the functional integrals can be expressed in terms of thermal 
averages of atomic observables. 

In section 5 first we analyze the effect of screening when X 3> A scree n by 
introducing the screened Coulomb potential. Various cases can occur according 
to the possible value of A scree n compared to the other lengths. When A scre en 3> 
Aph we show that the classical Lifshitz potential is exactly canceled by the 
thermal screening correction linear in T. Only diamagnetic terms beyond the 
dipolar approximation do survive in the leading tail of the atom- wall interaction 
(these terms will be presented in |20 j ). For densities such that either A screen <C 
Aph or even A screen <C A at , the same thermal screening correction linear in T 
has to be added to the Casimir-Polder or van der Waals potentials whenever 
Ascreen « ^ « A p h or A SC reen « ^ « A a t respectively. More comments are 
offered in the concluding remarks. 
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2 A generalization of the Feynmann-Kac formula 
to retarded interactions 



We consider an electron submitted to an external potential V(r) and in interac- 
tion with the quantum electromagnetic field. With q, p, [?' t ,p I/ ] = ifrSu U , /j,, v = 
1,2,3, the canonical quantum variables of the electron, the total Hamiltonian 
written in the Coulomb gauge is 

ff=^(p-^A(q)) 2 + V(q) + H iadl V-A = (4) 

The field is enclosed in a region A and obeys appropriate boundary conditions at 
the border of A. These boundary conditions, which do not need to be specified 
at this point, define a complete set of (real) orthogonal and divergence free field 
eigenmodes f 7 (r), J A drf 7 (r)fy(r) = if 7 ^ 7', V • f 7 (r) = 0. The eigenmode 
expansions of the vector potential and of the free photon Hamiltonian read 

A(r) = ^(a 7 + a 7 )f 7 (r) (5) 

7 

H ia d = ^2 e 7°7 a 7 ( 6 ) 

7 

where a 7 ,a 7 with [a 7 ,aL] = <5 77 ', are the photon creation and annihilation 
operators in the mode 7 with corresponding energy e 7 . Normalization factors 
entering into the definition of the vector potential are included in f 7 . 

The electron and the photons are supposed to be in thermal equilibrium at 
a common inverse temperature j3 = (/c^T) -1 (ks the Boltzmann constant), and 
we are interested in a functional integral representation of the effective electronic 
Gibbs weight when the field degrees of freedom have been traced out 

(e^rad = ^-IW-"* (7) 
^rad 

where Tr rac j denotes the partial trace on the Fock space of photon states and 
Zrad = Tr ra< ie~' 3ffrad is the partition function of the free photon field. To this 
aim, it is convenient to first introduce a functional integral representation of the 
action of the field variables. 

2.1 Path integral representation for photonic modes : the 
oscillator process 

For simplicity we deal first with a single mode f (r) of the field having energy e 
(dropping the mode index 7). We introduce the dimensionless canonical vari- 
ables (Q, P) of the mode, setting a + = ^2[3eQ, i(a^ — a) = [3eP so that 
[Q, P] — i. The dimensionless free Hamiltonian of this mode becomes 

(3H = (/3e)aU = \(P 2 + (f3e) 2 Q 2 - M (8) 
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and the contribution of this mode to the vector potential is A(r) = y/2{3eQ f (r). 

The trace over the states of one photon mode can be performed as an integral 
on the partial configurational oscillator matrix elements (R\e~^ H \R), where \R) 
is an eigenstate of the operator Q, so the partial average ([Tjl reads 

1 



dR(R\e~P H \R) 



(9) 



with Z = J dR(R\e~" °\R). To calculate these matrix elements we apply the 
method of path integral for bosons developed in chapter 5 of [T7], see also [26]. 
First, we split the total Hamiltonian as 

H = H el + H (10) 

where the electronic part 

H el = ^ (p - f VWeQ f(q)) 2 + V(q) = H A (Q) (11) 

depends both on the field operator Q and the electronic operators. The Trotter 
product formula applied to the partial matrix element (R\e~@ H \R) reads 

(R\e-P H \R) = (R\e-^-P H °\R) = lim {R\[ e -^ N e~^ N ] N \R) 



lim 

JV— >oo 



di?jv— l ■ 

f3H cl (R n )/N, 



dR n 



dRi e^ H ^' N (R\e'^' N \R N ^) ■ ■ ■ 

-m,i(Ri)/N( Rjp -f3H /N^ (12) 



xe —^" f {R n \e- ffH °/ N \R n _ 1 )---e- pu ^ til >^{R 1 \e 

A product of configurational matrix elements of the free oscillator 

(R\e-^'- SN ^ Ho \R N ) ■ ■ • (i?„|e- (s "- s '- l) ' 3Ho |i?„_i) • • ■ (i?i|e- (si - s)/3Ho |i?) 

(13) 

defines the oscillator process. The product ( [T3| is interpreted as the joint prob- 
ability (up to a normalization) for a closed path 9t(s) starting in R at time s to 
be found in R% at si . . . , in R n at s n . . . , and again in R at time s'. The ma- 
trix element (R'\e- ( - s ' - S ^ H °\R) is given by the Mehler formula (see [17] section 
1.8.1) which is a Gaussian function of R and R', so that the process defined 
by the distributions ( [l3| is GaussianJ^] The corresponding Gaussian functional 
integral is denoted by J^^Zr D[fK] 

Introducing the time dependent electronic Hamiltonian 



i 

2m 



WeX(s) f (q) + V(q) 



(14) 



associated with a given realization of the process, we see that the partial matrix 
element (12 1 can be represented by the functional integral on oscillator paths 



(R\t 



^ H \R) = 



m.(o)=R 



L>[5H]Texp 



-13 / dsH el (m(s)) 



lr The explicit form of Mehler formula will not be used here. 
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The need for the chronological ordering T in the propagator comes from the 
fact that H e \(d\(s)) is still an operator depending on the canonical variables of 
the electron. The partial average d9~} eventually reads 



Texp 



(3 / dsH el (m(s)) 



rad 



where 



)rad ^ s *° ^ e evama ted by the normalized integral 

r SK(l) = R 



' rad 



1 



dR 



D[m] 



(15) 



(16) 



SR(0)=fl 



when field quantities are expressed as functionals of the paths (Then Zq de- 

fined in (|qJ also reads Zq = j dR /Ju ).=ii D[9\].) The stationnary Gaussian pro- 
cess ( Jl6| is entirely defined by its covariance (!SH(s)!SH(s')) rad = — s 'l)^) ra d- 
The latter is easily calculated when we express it in operator form 



(X(sWs% ad = ±Tr [e-^T(Q(s)Q(s'))] 



(17) 



where Q(s) — e s ^ H °Qe S/3H ° is the imaginary time evolved operator. Using the 



commutation relations 

d 
ds 
d?_ 
ds 2 



Q(s) = [pH ,Q](s) = -iP(s) 
Q(«) = \fiH , {13H Q , Q]](s) = (/?e) 2 Q(s) 
one establishes that (9t(|s|)SH) rad obeys the differential equation 



ds 2 



<*(M)*>rad = <*(«) 



(18) 



(19) 



Its solution with periodic boundary condition (9t(s = 0)lH) rad = (5H(s = d 
is 

— /9e(X— , p -/3e\s\ 

<*(M)*>rad = T-/^ (2°) 



2/3e(l 



2.2 Path integral representation for the electronic variable 

We observe that the configurational matrix element of the electronic operator 
pl 



(r|{ e -^) rad |r) = ( (r|Texp 



f3 / dsH el {m(s)) 



|r) 



(21) 



rad 



is the field average of the imaginary time propagator associated with the Hamil- 
tonian H e \(?R(s)) (14 1 of an electron in presence of a classical time-dependent 
magnetic field, with potential vector A(r,s) = \/2j3eD\(s) f(r). We recall that 
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the Feynman-Kac-Ito representation of the configurational diagonal matrix ele- 
ment of this propagator , [T7] , [16] reads 



(r|Texp 



x exp 



-/3 / dsH el (m(s)) 
o 



(2ttA 2 ) 



2^/2 



-P / V(r(s)) 
Jo 



-P 1 ^ 



de(«) • f (r(a)) fR(s) (22) 



Here r(s) = r + A£(s) is a closed electronic path starting and ending in r. 
It is written in terms of a closed dimensionless Brownian path at the origin 

£(s), < s < 1, £(G) = f(l) = (a Brownian bridge) and A = is the 

thermal de Broglie length of the electron. The measure D[£] is the corresponding 
conditional Wiener measure normalized to 1. This measure is Gaussian, formally 
written as 



D[$ = cx p(-\f Q ds 



ds 



It has zero mean and covariance 



D[Z]em»(s') = S^(mm(s,s') - ss') , 



(23) 



(24) 



where £ M (s) are the Cartesian coordinates of £(s). In this representation a 
quantum point charge looks like a classical charged closed filament r(-) = (r, £) 
located at r and with a random shape £(s), Q < s < 1, the latter having a spatial 
extension given by the thermal de Broglie length (the typical size of quantum 
position fluctuations). The magnetic phase in ([22} is a stochastic line integral: 
it is the flux of the magnetic field across the closed filament. 



From the Feynman-Kac-Ito formula ( 22 I , the average over the radiation 



exp i 



degrees of freedom to be performed in (21 1 occurs in the form of the Fourier 
transform of the Gaussian measure whose covariance is (5H(s)K(s')) rad . Hence 
from the basic formula 

ds a(s) < R(s)]\ =cxp(-l- [ dsf ds 1 a(s) (9t(s)9t(s')) rad a(s') 

/ / rad V ^ Jq J q 

(25) 

for the Fourier transform of a Gaussian, the matrix element (21 1 eventually 
reads 

(r|( e -^) rad |r) = 

J D[£]exp (^fijds V{r(s))j exp (-^W rad [r,f]))(26) 



(2ttA 2 ) 
with 

Wrad [r,£] 



2A 2 ^ ^ 



de(s) dt»(s'){K(s)<R(s')) vad r(r( S ))r(r(s')) 
n c Jo Jo 



(27) 
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This formula can be readily extended to the case where the field has a multimode 
expansion ^ with mode variables Q 7 ,P 7 and mode energies e 7 . Reinstalling 
the mode indices and noting that the modes are independent and identically 
distributed, i.e. (9t Y (s)9 1 l 7 '(s')) ra d = <5 7 y (fH 7 (.s)5H 7 (s')) ra< i, Eq.(27l becomes 



o\2 /■! rl 

[ r ^=/^E J o <^( s ) j Q de(s')e^(^r(s)^(s')) Tad f!;(v( S ))f^T( S ')) 

(28) 

Eq. (26) is the desired generalization of the Feynman-Kac formula when the 
electron is not only submitted to a static potential V but also to the field-induced 
potential W ra d mediated by the presence of the photon field. 



3 Atom near a metallic wall 
3.1 Description of the system 

The system consists of a Hydrogen atom located in the vicinity of a metallic 
wall, the surface of which is located at x = 0, and interacting with a quantum 
electromagnetic field. The atom is made of an infinitely heavy proton at R = 
(X, 0,0), X > 0, and of a quantum electron of mass m and charge e at r = 
R + a, a being the relative position of the electron with respect to that of the 
proton. The wall is treated macroscopically : on the wall surface at x = 0, 
the electronic wave function is assumed to satisfy Dirichlet boundary conditions 
while the electromagnetic field satisfies the metallic boundary conditions. The 
electrostatic interaction with the wall is determined by the method of image 
charges. 

Thus the electrostatic potential part of the Hamiltonian |4} 

V(r) = - v ^-—-e 2 V im (r) (29) 
|r-R| 

is the sum of the Coulomb interaction between the electron and the proton in 
the bulk and of their interactions with image charges located at R* = (—X, 0, 0) 
for the proton and r* = R* + a*, a* = (—a x ,a y , a z ) for the electron: 

^) = \ (pnh*i + f^i - ~ F^ri ) (30) 

With the definition of the Fourier transform 

f(X) = J j^e-^fM (31) 

we can write 

1 + e ~ ik '( a - a *) - e ik ' a * - e - ik ' a l (32) 
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The large distance asymptotics of V5 m (X) is determined by the behaviour of its 
Fourier transform as k — * 0, which reads from (|32l 



Mm(k) 



2tt 



(k-a)(k-a*) 



(33) 



This implies that the atom and its image charge have the dipolar interaction 



e 2 V im (X) ~ - 



as X 



3(R-a)(R-a*) 



R=(2X,0,0) 



(2xy 



[4 + U a l + a D] 

(34) 



The electromagnetic field is enclosed in a metallic box L x L x L whose face 
at x — coincides with the wall position. This gives rise to the eigenmode 
expansion of the vector potential of the form ([5) with mode indices 7 = (k, rj) 



A"(r) 



k,?7 



L 3 



a:,y,z (35) 



where k are wave numbers and e^(k), 77 = 1,2, are two unit polarization vectors 
orthogonal to k. The functions /i^(r) are the cavity modes corresponding to a 
transverse electric field whose longitudinal part vanishes on the metallic faces 
of the box at x — 0, x = L, y = 0, y = L, z = and z = L: 

h x (k, r) = cos(k x x) sin(k y y) sin(k z z) 
h y (k, r) = sin^a;) cos(k y y) sin(fc z z) 
h z (k, r) = sm(k x x) cos(k y y) cos(k z z), 



h — u — 

L L 



y 



K = 



L 



n x ,n y ,n z = 0, 1,2, 



(36) 



In (351 g(A cut fc) is a smooth spherically symmetric ultraviolet cut-off, g(0) = 1, 
and g(A cu tfc) = when k > = mc/h. The prefactor in (35) includes the 
normalization of the eigenmodes ( |36| and assures the correct form of the free 
radiation Hamiltonian d6~} with eigenmode energies 



e fe = hu k , w k = ck 



(37) 



Note that in this model we have not included the spin of the electron and its 
Pauli coupling with the electromagnetic field. 



3.2 Effective atom-wall interaction 

Assuming that the electron and the photons are in thermal equilibrium at tem- 
perature T, we consider the immersion free energy of the atom in the photon 
field and in presence of the wall 

F(X,P)=-k B T lim m ,. - (38) 



Tr L je -" fl ™ 
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where the trace Tr L = Tr^TrJ^j runs over the space of electronic wave functions 
together with the states of the photonic Fock space which satisfy the boundary 
conditions on the box surfaces. The limit L — > oo means that the field region 



has been extended to the whole half-space x > 0. By virtue of (34 1 V\ m goes 
to zero when X goes to infinity ; henceforth, after the limit L — > oo has been 
taken, the bulk quantities can be obtained from the corresponding ones in the 
wall vicinity by taking the limit of an infinite X. Then the effective atom-wall 
interaction 

$(X,/3) =F(X,P)- lim F(X,(3) (39) 

X— >oo 

is defined as the difference between the immersion free energy F(X,f3) when 
the atom is located at distance X from the wall and its corresponding bulk 
value obtained as X — > oo. After inserting the definition ([7} into p8| and 
performing the electronic trace in the configuration representation, Q(X, (3) may 
be expressed as 

$(X,0) = -k B T (40) 



x In 



/ dr(r|(e-^( x )) rad |r)/ lim f rfr(r|( e -^W) rad |r) 

J X->+oo J 



It is understood here that the field extends in the half-space x > and the 
notation H(X) recalls that the Hamiltonian depends on the atom-wall distance 
X. Our main interest will be the analysis of Q(X, (3) as a function of the distance 
X for various temperature regimes. 

At this point we repeat the warning already given in the introduction about 
diverging atomic traces. The spatial integrals on electronic configurations in 



(40) must be restricted to some finite spatial region representing the effective 



available space for the atom in a low density phase. 

3.3 Functional expression of the effective atom-wall inter- 
action 



In view of (40 1 we have to specialize the generalized Feynman-Kac formula (26 1 
to the atom-wall model. The electrostatic potential takes the standard path 
integral form 

V(r,t) = [ dsV(r(s)) = -e 2 / ds 1 - e 2 / dsVi m (r(s)) (41) 
Jo Jo |r(s)-R| Jo 

whereas the effective magnetic potential VV^ ad (r,^) in the finite box L 3 becomes 

x g(A ph M-s')M k > r ( s ))M k > r ( s ')) (42) 
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This follows from a comparison between Q5J) and ( 35 1 showing that the mode 
functions in (28 1 (including the proper multiplicative factors) have to be iden- 
tified with 



/£>) 



IQ-kHc g{\ cut k) 



Vk 



<(k)ft"(r) 



In (421 



^„(k) = X;^CkK(k) = ^ 



(43) 



(44) 



is the transverse Kronecker symbol resulting from the summation on polarization 
indices and it is convenient to make the definition 



By virtue of ( 20 1 



Q(A ph M = (/^) 2 («k(M)<K k ) r ad 

A ph fc e-W(i-H) + e - WM 



(45) 



Q(A ph fc,s) 
A p h^ = P^k = flhck 



(46) 



The function Q(u, s) is the thermal photon propagator. As seen from (|46j it is 
the Green function of the quantum harmonic oscillator with periodic boundary 
conditions Q(u,0) — Q(u, 1). It satisfies the normalizations 



/ dsQ(u,s) = l, Iimfi(«,s) = l 
Jo 



(47) 



One has Q(u,s) > 0, Q(u, 1 — s) = Q(u,s), < s < 1 and Q(u, s) can be 
extended to a periodic function of s of period 1 for all s. 

It remains to extend the box L 3 to the whole half-space x > 0. As a result, 
lirnz,_, 00 W^ ad (r,^) = yV iei< i(x,£) will become independent of the location of the 
electronic filament £ in the y, z directions. So we can equivalently perform a 
spatial average L dy^J Q dz — , replacing sin (k y y) = 1/2, sin 2 (fc z z) — 
1/2, cos 2 (k y y) — 1/2, cos 2 (fc z z) = 1/2 when developing the products 
h^(k, r(s))h v (k, r(s')) in (|42|. We find for instance 



h x (k,r(s))h x (k,r(s')) 
1 



cos k x x(s) cos k x x(s') cos fc y [j/(s) — j/(s')] cos fc z [z(s) — z(s')] (48) 
Noting that 

cos k x x(s) cos fc^s') = § cos k x [2x + A£ x (s) + A^(s')] + \ cos k x \[£, x (s) - ^(s')} 

we see that the second term in the r.h.s is ie- independent. When L goes to infin- 
ity, the discrete sum on modes j-g ^ k where k x — irn x /L, k y — irn y /L and k z = 
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nn z /L with n x ,n y , n z =Q,l,2,... is replaced by integrals J +o ° ^f- J + °° J + °° 
Noting that <5^,(k) = fl — ^fj and the remainder of the integrand are even 

functions of k x , , k y and k z the integration can be extended over the whole of k 
space, e.g. 



dk 7 _ 



cos{k x A x ) 



dk 



and eventually 

lim - 

L->oo L 6 



[un ^ ^ £(^fi^ ( k ,r( S )) ^ (k,r(,s'))) 



(49) 



dk g(X cut k) 



(2tt) e 



e -ife J! [2a+A€ a! ( S )+Ae !c ( s ')] + g-ifc«A[^W-|«(t)] 



xe 



The second term in the r.h.s is ^-independent and thus contributes to the bulk 
value W^ d (£) = linx^oo W rad (x,£) of the field-induced potential. Dealing in a 
similar way with the other components, we eventually find 



W rad [x,£] = W r b ad [£] + W r w ad [ X) Q 



(50) 



where 



x e 



A 2 



2 r tik 



Q(A ph fc,s-s') 



(51) 



represents the bulk self-energy of the electron in the photon field. The expression 
of W™ d (a;,£) written in Fourier representation is 



W r w ad [z,£] 



dk 

(2^ 



-2ifc.-X 



(52) 



where J**" [a,k] is a functional of the position a(s) = a + A£(s) of the electron 
relative to the proton 

/^[a,k] = ^-/ da"(«) / da"{s')Q(\ ph k,s-s') 
A P h Jo Jo 

xg -ifc :c [Q x (s)+a x (s')] e -ife ! ,[Q !( (s)-Q ! ,(s')] e -ifc i [a i (s)-Q i (s')] 



and ( xx = 1, C M „ = — 1 otherwise. In (53 1 we have used the properties Ad£^(s) 



da„(s) and A[£ (s) - £(s')] = a(s) - a(s'). 



15 



Collecting these results in ( 26 1 , the explicit form of the generalized Feynman- 
Kac formula for the atom-wall system becomes 



(r|(e- 



) radlO 



2ttA 2 



3/2 



Be z / ds 
o 



0(X + a x (.))exp(pe 2 V im [X,a] 



2 



where 



^ m [X,a] 



dsV im (X,a(s)) 



(54) 



(55) 



The func- 



is the standard Feynman-Kac representation corresponding to (30 1 
tion 8(u) = 1 if u > 0, 6*(u) = if u < implements the Dirichlet boundary 
condition for the electron wave function at the wall. The second square bracket 
incorporates the effects of the wall and tends to 1 as X — > oo. Hence 

lim (r|(e-^W} rad |r) 

X— >oo 



( l 



3/2 



#(£)exp 



Be 2 I ds , } ci 



(r|(e 



/rad 



|r) 
(56) 



is the generalized Feynman-Kac formula for an atom immersed in an infinitely 
extended electromagnetic field in homogeneous space, with corresponding Hamil- 
tonian H^,. 

In fact the bulk self-energy Wi, d (£) of the electron in the homogeneous 
photon field is identical to that derived in [19] , formula (66), with the use 
of periodic boundary conditions. In other words, the generalized Feynman- 



Kac formula in the bulk Eq.(56l could of course have been established directly 



starting from the simpler Hamiltonian Hy, obtained by removing image charges 
and metallic boundary conditions at the wall right away. (In the method used 
in Ref.[19 for a quantum field in the absence of any metallic boundary, the 
photonic trace is expressed in the basis of the coherent states associated with 
the photonic modes, instead of the basis of photonic modes themselves, and 
then every corresponding matrix element of the Gibbs factor is replaced by a 
bosonic path integral different from that of the free oscillator process. However 
the result of the trace is the same, as it should.) The spectral properties of 
this Hamiltonian (called The Standard Model of Nonrelativist QED) have been 
extensively studied (see the recent paper |28] and references therein). H\> has an 
unique and non degenerate ground state, but atomic excited states are turned 
into resonances because of the phenomenon of spontaneous emission implying 
that the rest of the spectrum is continuous. 

Let us introduce the normalized atomic measure in the bulk associated with 
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JdajD[£]exp 






J da J D[(]exp 







(57) 



As discussed in the introduction, the J D^ t [a] integration needs a regularization 
by limiting the available space for the electron position in the da integral to a 
finite sphere. This regularization will be understood in D^ t [a] integrals without 



additional notation. Then, from (40 1, (54 1 and (571 the atom-wall effective 



potential receives its final form (for a fixed temperature T > 0) 
$(X) = -fc B Tln J Dl t [a] 9(X + a x (.)) exp (/?e 2 Mm [X,a] 



2 ' K rad 



[X,a]) 
(58) 



Since both T4 m [X, a] and W^ d [X, a) tend to zero as X — > oo, the large 
distance behaviour of the atom- wall potential is obtained by expanding ( 58 1 to 
first order in these potentials giving three contributions 



$(JQ ~ $ im (X) + $ rad (X) + $ geomc (X) 



<J> im (X) = -e 2 / D b At [a]V im [X,a} 



The two first ones 



^d(X) = -J C a >]W r w ad [XM 
arise from the image-charges and the photon field. The last one 



geomc 



(X) 



-k B Tln J Dl t [a]0{X + a x (.)) 



(59) 

(60) 
(61) 

(62) 



comes from the pure geometrical constraint imposed by the wall on the atomic 
weight. In view of the fact that Vi m [X, a] and W^j [X, a] vanish at large dis- 
tance, this constraint can be disregarded in (60 1 and (61 1 when looking for the 
dominant term as X — ► oo. In the low temperature regime considered in section 
4.3, where only ground state contributions of the Hydrogen atom will be kept, 
<& geomc (X) will be exponentially small for X ^> a B since the ground state is 
exponentially localized in the proton vicinity. 
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4 Single atom in a thermalized quantum electro- 
magnetic field 

4.1 The net interaction 

We examine the large distance behaviour of the Coulomb potential due to image 
charges ( 60 1 



$im(A) = 

and of the field-induced potential (|6T| 

*rad(*) = 



dk 



(2tt) 3 



dk 



$im(k) 



(2?r) s 



$rad(k) 



(63) 



(64) 



written in Fourier representation according to (31 1. As far as the Coulomb part 
is concerned, going to the path representation of (33 1 and averaging with the 
atomic weight D^ t , it follows immediately that 



$im(k) 



2ne 2 
W 
2Tre 2 



Dl t [a] / d s (k-a(s)(k-a*( s )) + 0(k 2 ) 



EW^ / D^[a]J ds [a x (s)] 2 + 0(k 2 ) 



(65) 



In the second line, we have used the rotational invariance of the bulk atomic 
measure to simplify the expression. Noting that ^2 (^^1 = 2k 2 — k 2 we 
conclude that the leading tail of the Coulomb interaction ( 63 1 arises from the 
singular part of $i m (k) which is 



$ im (k) ~ Aire 



k 2 
k 2 



D^[a) / ds [a x (s)Y 



(66) 



The small k behaviour of <J> rac j(k) requires a more elaborate study. According 
to (61 1 and (52 1 one has 

$rad(k) = ^f\j2 S U^ J D h Ja]I^[a,k] (67) 



with I^ v \a, k] defined in (53 1. The exponentials in 7 M "[a,k] can be expanded 
in powers of k 



k] = P lv [°1 [a, k] + J"" PI [a, k] 



(68) 



where the indices 0,2,... refer to exponentials expanded to 0,2,... order (odd 
orders will not contribute because averages of odd powers of a p vanish by rota- 
tional invariance of the atomic weight). Accordingly we have the expansion 



*ttd(k) = *^(k)+*N i (k) + 



(69) 
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This does not allow yet to infer the small k behaviour of $||^L(k) and $JSj(k) 
in a straightforward manner, since there is still a k dependence in the photon 
thermal function Q (X p hk, \ s — t\) and atomic averages have to be performed to 
obtain these quantities. The following facts turn out to be true: 

• Both ^l^jPO and ^JSLpQ contribute at the same leading order as X — > 
oo. 

• The term ^ad(^) corres P°nds to making the dipolar approximation, i.e. 
neglecting the phase occuring in the exponentials in (531. This approxi- 
mation is valid in an atomic state when the photon wave length is much 
larger than the Bohr radius, i.e. k ^> l/a B , or equivalently X 3> a B . 

• The next correction Q^JX) involves the diamagnetic susceptibility of 
the atom. Asymptotic contributions of ^^(X) are either \ C /X, with 
a B <C X <C A at , or af s smaller than those of ^^ A (X) [20], where ct{ s is the 
fine structure constant and A c is the Compton length A c = h/mc with 
Ac/as = <*fs- 



In the rest of this paper, we shall work in the dipolar approximation, re- 
stricting our attention to the first term <f>j.^ d (X), Consequently, from now on, 
our total potential will be restricted to the form 

$[°](X) = $ im (X)+$f a ] d (X) (70) 

[2] 

We shall present the analysis of the finer corrections arising from ^ T ^(X) m 
another paper |20j. 



Replacing the exponentials by 1 in ( 53 1 , we have 



I>* v M[a,k] = ^r [ da^(s) f da^(s')Q(X ph k,s- s') (71) 
\h Jo Jo 



where Q(A p hfc, s — s') is in fact a function of \s — s'\. Rotational invariance of 
the measure D^[a] implies that / D% t [a] 1^ [a, k] M = S^ v I{k) with 



ph 



( Dl t [a] f da x (s) [ da x (s')Q(X ph k,s-s') (72) 
J Jo Jo 



where we have used the symmetry of the integrand under the exchange of s, s' 
to reduce the time integration to the sector s' < s. From the definitions (61 1, 
(64) and |52|), and in view of the relation £ w (l - k^/k 2 ) = -2kl/k 2 , we 
get eventually 

^) = -^f 2 I f (73) 
All the relevant information is now contained in the function I(k). 
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Then one can get rid of the stochastic integrals by means of a double in- 
tegration by parts. The calculation can be formally performed by using the 
standard rules for ordinary differentials. (A more detailed justification for this 
calculation will be found in [20 ). When using standard integration by parts, 
one takes into account the fact that for a Brownian brigde a(l) = a(0) and the 
following properties of the function /(it) = Q(X p hk,u) arising from the defini- 



tion (461 : /(l -u) = /(it) and /'(l - u) = -f'(u) for < u < 1. As a result, 



the expression d72| is equal to 



Ay(fc) = -2 J Dl t [a]Q'(\ pil k,Q) j\s [a x (s)} 2 (74) 

-2 f Dl t [a] f ds f dtQ"(X ph k,s-t)a x (s)a x (t) 
J Jo Jo 



The expression (46 1 of Q{k,u) leads to 



Q'(k, 0) = -^Ap h /c 2 and Q"(\ p hk,u) = Xp h k 2 Q(X ph k,u), u>0 (75) 
so that 

I -^- = fds / D atW M*)f - (76) 

with 

B(fc) = 2 / ds / dtQ(\ ph k,s-t) [ D b Ja}a x (s)a x (t) (77) 
Jo Jo J 

Finally inserting the expression I(k) into ( |73| | gives 

if a ] d (k) = -4^ 2 M I' d S J D b at [a] [a x (s)} 2 + ^e^B(k) (78) 



At this point one observes the remarkable fact that the first term in (78 1 is 
nothing else, up to the sign, than the singular small-k expression ( |65| of the 
Coulombic part of the interaction. Consequently, when considering now the 



total potential ( 70 1 the dipolar long- ranged part ~ X 3 due to <j>; m (X) is exactly 



compensated. Thus, the total net potential at large distance 

~ 4.e 2 J 'j^e^ f 2 B{k) (79) 



is entirely determined by the behaviour of the function B{k) (77). Clearly this 
function shows an interplay between the photon thermal propagator and the 
squared moduli of the atomic moments, which will be made explicit in the next 
section. 

4.2 Asymptotic regimes 

Although we have disregarded diamagnetic contributions, there are still effects 



depending on the fine structure constant in the expression (79) of <5>( )(AT). In- 



deed since our calculation so far is not perturbative with respect to the coupling 
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constant (the charge of the electron), the thermal weight -D at [c*] (57 1 of the atom 
still includes its coupling with the electromagnetic field. When this coupling is 
switched on, it turns the excited states of the bare hydrogen atom into reso- 
nances, with displaced energy levels and finite life times due to the phenomenon 
of spontaneous emission (see comments after (56l). 

We assume that these effects do not modify the power of the X-decay but 
only bring tiny corrections to its amplitude. Therfore we neglect them in the 
subsequent analysis of ^°~>(X). This amounts to drop the effective interaction 
W^ ad [£] in (571, namely replacing Dg t [a] by the bare hydrogen atom weight 



Coming back to the operator formulation 



/da/D[£]exp 


\_(3e 2 J Q rfs^_ 




J da J £>[£]exp 


f/3e 2 / dS| Q( 1 s) |_ 





(A) at = -^TrU 
^at 



-/3-f/a 



A) 



Z»t = Tr 



H a t — 



P 

1m 



|r-R| 



(80) 



(81) 



is the usual thermal average for an atomic observable A. It is important to 
repeat here the warning about our phenomenological treatment of screening 
(see the introduction and the remark after (57 1). The above formulae have to 
be regularized, either by a cut-off in the spatial da integral occuring in (80 1 or 
by an energy cut-off limiting the evaluation of the traces (81 1 to a finite number 
of hydrogen eigenstates with maximal energy E max . 

With this limitation we can proceed to an explicit calculation of the function 
B(k) (77 1, which now reads 

B(k) = 2 [ ds f dt( a;c (s)a x (t)) at Q(A ph fc,s-t) 



(82) 



Coming back to operator langage by inverse Feynman-Kac transformation, the 
atomic fluctuation can be evaluated in a basis of eigenf unctions for the Hamil- 
tonian iJ at 



{a x (s)a x (t)) &t = 



1 



-Tr 



EE 



l(*Kb')l ; 



-[l-(B-t)]0Ei -{s-^pEj 



(83) 



where i max is the uppest state index such that Ei is equal to the phenomological 
cut-off E m3x discussed in the introduction. The rotational invariance of iJ at 
enforces that (i|a x |i) = so that the term j — i can be omitted [^J According 
to the expression p6| of Q(X p h_k, s — i), we have to calculate 



J a.t 



(3 / ds 



n 



dt 



(i 



-(s-t)/3e k 



-[l-{s-t)]0e k 



-0Ei .-{a-tWEj-Ei) 



(84) 



2 The j-summation on intermediate states runs over the whole spectrum of tha Hydrogen 
atom including its countinuous part 
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with e k = hck. The result of the integration is the sum of two terms 7jj + Jjj 



J A 



1 

-^at 
1 



£y — Ei 



£y — Ei — e k 



and therefore 



e -/3(Bj+e fc ) _ e -/3B, e -/3(£;+£ fc ) 



(£,■ - + e k y 



-0Ej 



{E 3 - Ei - e k y 



-Jji 



B{k) 



Eft 



(1 - e-^) Z, 



i\a x \j)\\h 



(85) 
(86) 

(87) 



i j^i 

We discuss now various distance dependences in the low-temperature regime. 
4.3 Low-temperature regime 

We specialize now our study to the low temperature regime (JT|. In this situ- 
ation we neglect exponentially decaying contributions of the order e -0( E i- E o) ; 
which amounts to only keep in (87) the terms having at least one ground state 
contribution [29] [22], i.e Z at = e^ E ° (l + 0(e-^ El ~ B ^) and 



B(k) = 



/3E 



(1 - e-^) 



]T \(0\a x \j)\ 2 (I oj + J 0j ) + ]T |(z| a;c |0)| 2 (/, + J l0 ) 

j?0 i 



(88) 



As explained in the introduction, once only the contributions with at least one 
state equal to the ground state are retained, the double sums becomes convergent 
and we can remove the cut off i max in the second sum of (88 1. Then sums 
involving the J function compensate each other because of the antisymmetry of 
Jij (see (86 1, and also |(0|aj;|j)| = |(j|aa:|0)|, and we are eventually left with 

B LT (k) = ]T |(0|a,|i)| 2 G(^ - E ,e k ,(3) (89) 
where G(Ej — E ,e k ,f3) is the following function of energies and temperature 



G(E j -E ,e k ,/3) 



Eft 



(1 



Ej — E + e k Ej - E Q - e k 



(90) 



The superscript LT in (89 1 denotes the low-temperature regime where 
exp [—(Ex — E )/k B T] terms are disregarded. 

Subsequently, in the X-behaviour of <J>[°1 (X) in the latter temperature regime, 
it is consistent to neglect the oscillating function of 2X/X !it with a damping fac- 
tor exp [— (Ei — Eo)/k B T] which arises in the inverse Fourier transform (79 1 
from the singular term e~^ ek /(Ej — E — e k + i x + ) (where Ej — E > and 
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e k > 0.). We recall that A a t = hc/(Eo — Ex) is the wavelength of the photon 
emitted when the atom jumps from the first excited state. 

To distinguish now the various possible large distance tails in the low- 
temperature limit it is useful to make the scaling k = q/X in the Fourier integral 
(79) where q is a dimensionless Fourier variable, leading to B(k) = B(q/X) and 

(91) 



e k = € q/X. Hence in (90 1 we have 



and 



Ej — E ± e k — Ej — E ± 



X 



hcq he f X 



X 



X V A. 



±q 



0e k = 



(92) 



Since A a t/A p h = ksT/(Ei — Eq) <C 1 in the considered low-temperature regime 
(JlJ, one has to discuss separately the cases when X is much smaller (much 
larger) than A a t or A p h- 

4.3.1 Electrostatic dipole interaction X <C A at <C A p h 



According to (91 1, when X <§; A a t we can neglect Ej — Eq compared to e k in 
d90|, so that 



G{Ej-E ,e k ,(3) 



1 + 



1 



KKtk \X at k / 
Subsequently, since £\ l(°l a x|j)| 2 = (0ja,|j0), Q simply yields 

B^{k) = (0|<£|0) 

l<gA at fc 



By the inverse Fourier transform (|79 1 



1 1 

A at «x X 3 4 



e 2 (0| a 2|0) = ^ dW (JC) 



(93) 



(94) 



(95) 



Thus at distances X <C A at we recover the standard electrostatic dipolar inter- 
action, which could have been inferred from (34) by replacing there the a 2 by 
mean ground state atomic moments (0|a^|0). 

4.3.2 Retarded interaction A a t <C X <C A p h 

When A a t <C X we can neglect e k compared to Ej — Eq in ( 90 1 



G(E j -E ,e k ,P) = S Eo ,E j +\^ - 5 Ei ,E 3 ] efccoth 



2 



1 

Ej — Eq 



[l + G(A at fc)] 
(96) 



Moreover for X <C A p h the argument (3e k /2 of coth can be considered to be very 
large (see (92 1 ) 



e fc coth (X ph k) = Hck [1 + (e- A p" fe )l 



(97) 
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and inserting (97 1 and (96 1 in (89 1 



v-I(0KIj)I 2 



B™(k) = hck^ y ^ ee B CP (k) (98) 

A at fc<l f—* £j; — Hin 



(The ground state is nondegenerate so that the latter sum involves no singular- 
ity.) The corresponding spatial decay follows from the inverse Fourier transform 



(791 



1 j A at «x«A ph IMtt^ Ei-Eo X*8tt aE ~ CP[ ' iyyj 



where 



o 2 y l(0Klj)l 2 , inn . 



denotes the static polarizability of the hydrogen atom in its ground state (the 
response of the atom to a classical external electric field in a given direction). 
This is precisely the formula originally found by Casimir and Polder [7] for 
retarded interaction in the atom ground state. Our derivation does not involve 
explicitly the concept of retardation associated with the propagation of Maxwell 
waves as in [7]. The change in the decay regime from X~ 3 to X~ 4 occurs as a 
consequence of the behaviour of the photon thermal propagator Q(X p hk,t) for 



large wave number k producing the extra k factor in ( 97 1 and ( 98 1 . Note that 
the original ground state Casimir-Polder formula remains valid as long as one 
neglects thermal effects 0(e~^ El ~ B °^). 

4.3.3 Classical field regime A p h -C X 

When X ^> A p h, /3efe -C 1 and we expand coth for small argument 

e fe coth(A ph fc) = 2k B T [1 + O (A ph fc)] (101) 



with the result 



B-(fc) = 2 fc B T^ fetf ee B class (k) (102) 



and 



$[°]m ~ -^ fi TV I(°I°»IJ)I 2 = -^-k B T a E = <& class (X) (103) 
V ' A pil «x 2Jf 3 B ^ Ej-E 4X 3 S B classV ; V ' 

The same result can be obtained by expanding the photon propagator for small 
k. Since Q(\ p hk,t) is an even function of k it behaves as Q(X p hk,t) = l + 0(k 2 ) 



(see (46 1). The term of order k in B(k) (|82j gives an analytic contribution to 
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VdW 



CP 



X 



Figure 1: Low-temperature atom-wall potential when A scree n is larger than all 
length scales at stake. 



Q(A ph M) = 1 in (821 



the integrand in (79 1 and rapidly decaying terms as X — > oo. Thus we can set 



-Bciass(fc) = 2 / ds dt(a x (s)a x (t)) at 
Jo Jo 



(104) 



or equivalently set k = or A p h = in subsequent formulae, which leads to 
(1031. Setting A p h = (3hc — is the same as treating the field classically by 
turning off the Planck constant in field expressions, hence (1031 describes the 



atom-wall force when the atom is immersed in a classical electromagnetic field. 
This potential vanishes at T = 0. 

We stress again the subtle behaviour of the force at finite (but low) temper- 
ature: three successive ranges of decays occur, as shown in Fig.l, 



X~ 6 if X < A a 



x- 



if A at < X < A 



P h. 



X~ 6 if X > A 



P h 



Note that the asymptotic formulae (951, (991 and (1031 are exact in the sense 



that they do not depend on regularization procedures of Coulombic traces. This 



is because in (881 one has only retained contributions of the ground state, which 
is localized in space. 



5 Screening effects 

As already said, as soon as the temperature is different from zero, the concept of 
isolated atom does not make sense. One must rather consider a non zero density 
equilibrium phase of nuclei and electrons: if both the density and the tempera- 
ture are sufficiently low, the latter can predominantly be found in atomic bound 
states, but there always remains a fraction of dissociated charges that provide 
a screening of the bare Coulomb interaction on distances greater than A scree n- 
This situation is precisely described by the so called atomic limit that defines 
the Saha regime (See Ref.|25] (section 7) and Ref.[22j). One obtains another 
case when the atoms under consideration are also in equilibrium with other 
species of dissociated charges, like in a strongly ionized gas or an electrolyte. 
The simplest model of the latter case is obtained by embedding the atom in a 
classical weakly coupled plasma. Here again Coulomb interactions are screened 
beyond some typical length A scree n depending on the plasma density. In both 
cases, screening has two main effects on the atomic phase 

• Regularization of divergent atomic traces 
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• Reduction of the range of inter-atomic forces 

Concerning the first point, we just recall here that the regularization consists 
in substracting to the Coulomb weight in ( |80| a number of terms of its large 

distance expansion, namely replacing exp fie 2 f Q ds i^feyi by 



exp 



1 



fie 2 / ds^ 



l«0)l 



E 1 

^— ' n 

n=0 



\ot(a)\ 



a(s)=a + \£(s) (105) 



which is integrable at a = oo whenever N > 3. The equivalent truncation in 
operator language consists in substracting the first terms of the Dyson expan- 
sion of the thermal propagator (see (5.11) and (5.12) in [22] or (117) and (119) 
in |30j). This truncation is by no means arbitrary, it follows from the system- 
atic treatment of screening provided by the so called screened cluster expansion 
presented in [21] . Once these truncations have been introduced, it can be rig- 
orously established, following section 5.2 of [22] or section 6.2 of [30] and the 
appendices of these papers, that the remainder in (881 is indeed 0(e^ /3 ^ El ^ E °' > ) 
(up to a polynomial in (3). 

Concerning the second point, the bare Coulomb potential between two charges 
ei , e 2 (written in Fourier representation) 

^c(k,6,6) = f dse^^W-A^M) (106) 

fc Jo 

with corresponding paths ri + Ai£i(s), r2 + A2^2(s) becomes a screened effective 
potential 

V SCIeen (k,fi, 6) = ^en(k,£l,&) + ^cr S ee„(k,£l,6) (107) 

which is the sum of two contributions [24] [22] (section 3.2). The first one 

fe„(k,6,6) = jf'd^c^w-^w) (108) 

has the Debye-Hiickel form familiar in the classical theory of screening. The 
replacement An/k 2 — ► 47r/(fc 2 + k 2 ), with k = A^J^ the inverse screening 
length, leads to an exponentially fast decay on scale A screen as |r x — r 2 | — > oo. 
In the limit of high atomic dilution, the second contribution takes the form 



[ d Sl [ ds 2 [5( Sl - s 2 ) - 1] (k • Aifi(«i))(k • A 2 6(s 2 )), k - 
Jo Jo 



47reie2 



(109) 



The spatial decay ~ |ri — r 2 \ 3 of the corresponding potential is dipolar and 
has a pure quantum origin. It represents the interaction of the two fluctuating 
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dipoles Ai£i(si) and A2^2(s2) generated by the intrinsic fluctuations of quantum 
charge positions. Here the screening is non exponential, causing only a reduction 
of the bare Coulomb decay <~ |ri — i~2 1 1 to the dipolar one <~ |ri — r2 1 — 3 . 
This algebraically decaying term disappears whenever one or both charges are 
classical. A more thorough discussion of quantum screening effects can be found 
in |31l 1321 1251 133] . When |ri — 1"2 1 3> A S creen one can disregard the exponentially 
small contribution of V^reen but one must of course keep the long range part 

T?alg 
v screen- 

In our model, when X ^> A screen , this entails replacing the bare Co ulomb 
interaction of the electron with its image by V^ r | en and thus, according to ( 109 1, 



replacing ( 65 I by 



<J> al s fkl = 

screen v'V 

— J D b Ja] j ds x j\s 2 [S( Sl -s 2 )-l] (k-a( Sl )(k-a*(s 2 ))+0(k 2 ) 
= $ to (k)-^f / Dl t [a] J d Sl J ds 2 (k-a( Sl )(k-a*(s 2 )) + 0(k 2 ) 

(HO) 

Using rotational invariance and neglecting fine structure constant effects in the 
atomic measure as before, one recognizes that the second term is identical to 



-2Tre 2 (kl/k 2 )B class (k) (see ([104} and ([102]). Thus the net potential ([79} is 
modified to 

f ,71, 1.2 

$M een (X) ~ ine 2 J —e~^ 2X f 2 \B(k) B clasa (fc)] = ^\X) - $ class (X) 

(HI) 

when X ^> A scree n and remains unchanged when X <C A scree n- We summarize 
below the long distance behaviour of the atom-wall potential according to the 
value of A scre en compared to the other lengths A pn and A a t- 

• A a t <C A p h <C A scr een- If X <C A SC reen the behaviour of is that given 



in su bsec tions 4.3.1, 4.3.2 and 4.1.3. If X A SC reen, we see from (111 I 



and (1031 that there is an exact compensation , &°'(X) ~ 0. Charges and 



thermalized photons conspire to cancel the leading order term ~ X 



-3 



A a t <C Agcreen "C A pn . If X -C Agcreen, &°'(X) behaves as in subsections 
4.3.1 and 4.3.2. However one has 

~$ CP (X)-$ claS9 (X), Ascreen « X « A ph (112) 

and - when X > A ph . 
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^ I ' I ' I ' ^ 

^at ^ph ^screen 
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(b) 
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I I I — 

■^screen ^at -^ph 
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Figure 2: Low-temperature atom- wall potential when A screen is of the same order 
as some characteristic length scale of the microscopic model. 



• A screen <C A at <C A p h- If X <C A screen , <&[°1(X) behaves as in subsection 
4.3.1, but 

screen 

$[°](X)~$ CP (X)-<& class pO, A at «X«A ph (113) 
and - when X > A ph . 

The various formulae are summarized in Fig. 2a Fig.2b and Fig. 2c. where <J>^ 2 ^ 
denotes the diamagnetic contribution to the decay mentioned in subsection 4.1. 

are analogous 
They 



The thermal corrections in ( 1 12 I and (1131, linear in T (see ( 103 1 ) 
to those found for van der Waals potentials at finite temperature in [22 
account for free charge screening, and because of the negative sign, leading to 



a weakening of the standard dipole (95 I and Casimir-Polder (99 1 interactions. 



6 Concluding remarks 

In this paper we have developed a functional integral method to analyze dis- 
persive forces at a microscopic level taking retardation effects into account. 
This formalism enables one to extract the asymptotic form of the force with- 
out recourse to perturbation with respect to the electromagnetic field coupling. 
Miyao and Spohn j3l] have applied the same method to give a non-perturbative 
derivation of van der Waals forces in the atomic ground states. Here we have 
considered the atom-wall forces at non zero (but low) temperature and shown 
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in a first stage how to recover van der Waals, Casimir-Polder and classical Lif- 
shitz forces in the dipole approximation. We have pointed out a consequence of 
screening that is not obtained in the usual applications of the Lifshitz theory. If 
the Coulomb potential between the atom and the wall mirror charges is screened 
by ionized electrons or possibly other types of mobile charges, we find a tempera- 
ture correction to the van der Waals or Casimir-Polder force which is linear in T. 
This correction originates from the fact that the screening of quantum charges 
is not exponential, but algebraic, and thus participates in building up the long 
tail of dispersive forces. Algebraic screening entails non analytic terms in the 
small wavenumber expansion of the dielectric function e(k,cj) [35 . Therefore 
applications of the Lifshitz theory that use simple analytic forms of dielectric 
functions to describe the different media do not predict the modifications of 
dispersive forces due to algebraic screening. 

Moreover the spontaneous emission of a photon by an excited state with 
energy Ei gives rise to cos[2Xhc/(Ei — E )] oscillations with a damping 1/X 
factor in the large-distance behavior of the atom-wall interaction. Such a decay 
comes up in our treatment (see the remark after ( 90 1 ) as well as in the study 
of an atom prepared in an excited state f5] or in the calculations for a sim- 
plified model for an atom with only two energy levels [12 , as already pointed 
out in the introduction. This tail, which has an exponential thermal weight 
cxp [—(Ex — E )/k B T], has been disregarded in the low-temperature limit con- 
sidered in the present paper. However spontaneous emission effects should be 
retained in thermal corrections whenever such exponential contributions can no 
more be neglected at higher temperature. 
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